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Abstract 



The anomalous behavior of the spin-spin correlation at large momentum 
transfer in pp elastic scattering is described in terms of the degrees of freedom 
associated with the onset of the charm threshold. A non-perturbative analy- 
sis based on the symmetries of QCD is used to extract the relevant dynamics 
of the charmonium-proton interaction. The enhancement to pp amplitudes 
and their phase follow from analyticity and unitarity, giving a plausible ex- 
planation of the spin anomaly. The interaction between cc and light quarks in 
nuclei may form a distinct kind of nuclear matter, nuclear bound quarkonium. 
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The fundamental description of the strong forces in terms of structureless constituents 
which carry the symmetries within the hadrons, and their interaction in terms of gluonic 
fields, is the basis of quantum chromodynamics, QCD. Despite considerable success in the 
application of QCD in the perturbative regime (PQCD) [|l|, there are important discrepancies 
with measurements related to the production and decay of heavy quark systems and 
the study of spin effects at high energy f^. A notorious example is the unexpected result 
observed in the polarization correlation Ann found in proton-proton collisions at large energy 
and momentum transfer ||^ . At ^/s near 5 GeV and 6'c.m. = 90° the rate of protons scattering 
with incident spins parallel and normal to the scattering plane, relative to the protons 
scattering with antiparallel spins, rises dramatically with Ann reaching a value of 60 % in 
a region where scaling behavior is expected. 

A simple explanation for this anomaly has its origin in the observation that at ^/s = 
5 GeV new degrees of freedom associated with the onset of a heavy quark-antiquark thresh- 
old could invalidate the scale-invariant nature of the basic constituent interaction at high 
energies by the interaction of light and heavy quarks at low relative velocity f^. Other mod- 
els had also been advanced to explain the spin anomaly, but they imply a radical departure 
from the PQCD framework or the Standard Model 0, introducing extraneous dynamical 
effects. The phenomenological analysis of Ref . P| , based onJ = L = S' = l broad resonance 
structures at the strange and charm thresholds interfering with PQCD quark-interchange 
amplitudes, successfully describes the anomalous p\ah dependence of the spin correlation but 
lacks of a clear dynamical content since the parameters of the model are chosen to fit the 
data. 

It is one of the main objectives of the present paper to show that the enhancement of 
the proton-proton elastic amplitudes close to the charm threshold is obtained from precise 
knowledge of the dynamics of heavy quarkonium with light hadrons. The spin correlation 
close to the threshold is determined from unitarity and analyticity without introducing 
arbitrary parameters and is largely predictable from basic principles. 

The study of the interaction of quarkonia with light hadrons is of particular relevance 
to a better understanding of the dynamics of gluons from a new perspective. Since dif- 
ferent quark flavors are involved in the quarkonium-nucleon interaction, there is no quark 
exchange to first order in elastic scattering and multiple gluon exchange is the dominant 
contribution p[. Remarkably, this system involves basic QCD degrees of freedom at the 
nuclear level, in contrast with the usual description of the dynamics of the nuclear forces 
with phenomenological lagrangians in terms of meson and nucleon degrees of freedom. 

Due to the small size of charmonium at hadronic scale, its interaction with external gluon 
fields from interacting hadrons is expressed as a QCD multipole expansion with higher 
order terms suppressed by powers of the quarkonium radius over the gluon wavelength, 
resulting in a long distance effective theory that has been applied successfully at the leading 
order in many processes involving low energy gluons [|T0|. The leading term in the expansion 
is the bilinear chromo-electric term, which is in fact the leading local gauge-invariant operator 
in an operator product expansion (OPE), whose matrix elements between hadronic states 
are evaluated using the trace anomaly ||lT],|12[. This is a notorious example of a calculation 
where QCD matrix elements, otherwise untractable by perturbation theory, are determined 
by the symmetries of the theory. 
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Using non-relativistic normalization of states and neglecting the charmonium recoil, the 
charmonium-nucleon amplitude can be written as a product of two matrix elements, one 
depending on the internal charmonium degrees of freedom and the other containing the 
external gluon fields 0. After some color algebra the result is 

where (H^^^ + e)~^ propagates the cc system in the octet state, uij-ed is the reduced mass of 
the proton-charmonium system, K its relative momentum and ip is the internal charmonium 
wave function. Standard evaluation of the charmonium matrix element gives 

'^red _7 3/1 



'^2ao(oEla-Ea), 



271 "'2 

where d2 = 287r/27 is the Wilson coefficient for a Is charmonium state with Bohr radius Oq. 

The field operator is expressed as a linear combination of the gluonic contribution to 
the energy-momentum tensor T^i'^^^^^ and the trace of the total energy-momentum tensor, 
T^, through the scale anomaly |[TT|Jl^ renormalized at the charmonium scale Aq ~ l/ao- 



At zero momentum transfer the matrix elements of T^i^^j^ic ^'^d are expressed in terms 
of the gluon-momentum fraction in the nucleon, V2, and the nucleon mass, M, respectively. 
In this case the value of / equals the scattering length, a^, and is given by 

Using in the previous equation a value of V2 = 0.5 from deep-inelastic scattering measure- 
ments, the hydrogenoid charmonium value ao = 1.6 GeV~^, the value of as = 0.6 evaluated 
at the scale Aq, and approximating the Gell-Mann-Low beta function by its leading term 
— 9a^/27r, we find a value of = —0.2 fm for the charmonium-proton scattering length, 
which agrees with the result obtained by Brodsky and Miller from Ref. [0. The sign of 



the amplitude corresponds to an attractive interaction. Other results [|T3|-[T^ differ widely 
because of different normalizations used in the calculation 

Scattering at low energies is described in terms of two parameters, the scattering length 
tts, and the effective range r^. We use the value of the scattering length determined above, 
ttg = —0.2 fm, to describe the low energy charmonium-nucleon scattering. The value of the 
cc-p effective range is unknown and can not be obtained from the previous calculation, valid 
only at zero momentum transfer. We expect however a value for re between ao — 0.3 fm 
and 1 fm, which accounts for the small size of charmonium and the short range of the QCD 
van der Waals force. With these values, we compute the effect on the pp elastic amplitudes 
at the onset of charm threshold and determine the binding energies of charmonium in 
nuclear matter |]7|JT7[]. Our results do not depend significantly on the value of rg. 

The unitarity condition relates the imaginary part of an elastic scattering amplitude to 
a sum over all possible intermediate states that can be reached from the initial state. The 
enhancement to a given partial-wave amplitude T is determined by the analytic properties 
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of the transition amplitude in the u = complex plane, with p the cms momentum, and is 
given by the dispersion integral 



^ ^ ^ ^ -kJq v'-v-ie 

We assume that the inelastic contribution of unitarity to the pp elastic scattering amplitude 
is dominated near threshold by the creation of a ppcc state with low relative momenta and 
internal quantum numbers j = 1, / = 0, which corresponds to the J = L = S = 1 quantum 
numbers of the pp system. The triplet wave has maximal spin correlation Aj^^ = 1, and 
the pp orbital angular momentum must be odd since the cc has negative intrinsic parity. 
The ppcc state is produced with partial width characterized by the measured branching 
ratio = 0.0155. The enhancement to the partial wave amplitude Tjj = Tj^l can be 
represented in the N/D form in terms of the discontinuity of T(z/) 



^•^ ^ ^ D{u) 2711 J v-v< 
where the Jost function, D{y), is defined by 

1 r b{y')dv' 



D{y) = exp 



TT J z/ — z/' — 



in terms of the proton-proton-charmonium scattering s-wave phase shift (5(z/) which, by 
unitarity, is the phase of the pp enhanced amplitude. We limit our calculation to the 
zeroth iteration D = 1 in the previous integral, to obtain the determinantal approximation 
rpenh _ J" I D = JIT with R = D^^ the enhancement factor. 

The pp elastic scattering is conveniently described in terms of the usual five helicity 

amplitudes: 0i = 02 = ( 03 = 04 = (-+I0I+-) and 

05 = (++|0|H — ). At large momentum transfer the PQCD quark interchange model (QIM) 
amplitudes |20| for proton-proton scattering dominate over quark annihilation and gluon 



exchange [^. The helicity conserving amplitudes are related by 0^^^ = 20^^^ = —20^^^ 
and are the only non- vanishing at leading order We assume the same form for the 

QIM amplitudes as in Ref. expressed in terms of the standard proton dipole form factor, 
consistent with nominal scaling law and angular distribution We do not consider any 
phase dependence on the QIM amplitudes since they cancel out in the present calculation. 

Since the nominal power-law scaling for the scattering amplitude is not preserved indi- 
vidually by each partial wave, a modification of a given partial amplitude effectively breaks 
the s""^ scaling behavior of the total amplitude at the charm threshold. The partial wave 
expansion of the helicity amplitudes is expressed in terms of singlet, coupled and uncou- 
pled triplet amplitudes. Retaining the enhancement in the J = L = S = 1 partial wave 
amplitude we need only to modify the J = 1 wave in the uncoupled-triplet partial wave 
expansion 

03- E {2J+l)Tjjdi, 

J=odd 

04- - E {2J + l)Tjjdi,,, 

J=odd 
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with the Wigner coefficients. The enhanced hehcity amphtudes (jf"^^ expressed in terms 
of the QIM hehcity amphtudes cf)^™, the branching ratio i^cf, and the s-wave enhancement 
factor R{v) are 



/^enh ^ ^QIM ^ ^Bl?{R - l)TnrfJi, 
^enh _ ^QIM _ ^BIP{R - l)Tnrf^i, 



where Tn = ^ /-^i^('iii03 — '^iii04 and 02°*^ and 05°*^ are zero to leading order. 

For the pp-cc state the s-wave phase shift is expressed in terms of the effective range 
approximation as kcotS = —l/a + ^rk"^, where k = p — po, with po = 2.32 GeV/c the 
cms momentum corresponding to the charm threshold. In this approximation the Jost 



function is obtained directly from the phase shift |18| giving for the enhancement factor 
R{k) = {k + ia)/{k + where a and (3 are related to the scattering parameters by 
a = (a — [3)/aj3 and r = 2/(a — [3). For definiteness we take a = 2as, r ~ r^, assuming 
simple additivity in the interaction of charmonium with the external fields of the two protons. 
Comparison of the prediction for the spin correlation A^j^ at 9c.m. = 90° with the available 
data near the charm threshold, is shown in Fig. |l] for different values of rg. 

Since the QCD van der Waals force is attractive []7|,p!2|- |l5|] , it could lead to the formation 
of quarkonium bound to nuclei [[F|,p!7|. The discovery of such state would unveil the purely 
gluonic component of the strong forces, a genuine QCD effect at the nuclear level, which do 
not involve meson degrees of freedom. In fact, in the absence of valence light-quark exchange 
there is no one-meson standard nuclear potential and the contribution of higher-order inter- 
mediate meson states is negligible Furthermore, there is no short-range nuclear repulsion 
from Pauli blocking. In Ref. an estimate of the QCD nucleon- charmonium van der Waals 
potential was obtained by rescaling the high-energy meson-nucleon interaction described 
by the Pomeron model [|^, and determining the coupling and range of the interaction 
parametrized by a Yukawa form, V{r) = —'ye~^'^/r, by adjusting the meson form factor in 
the Pomeron amplitude to describe the size of charmonium. The scattering length in the 
Born approximation, 2'ymj-ed/ fJ''^ — 0.5 fm, overestimates the strength of the interaction by 
almost a factor of three compared to the QCD non-perturbative result obtained above. 

Due to the short range of the nuclear forces, additivity is not valid for all the nuclei but 
the very light and the interaction with charmonium depends on the nucleon distribution. 



Following Wasson |T^, we write the charmonium- nucleus potential, Vcc-Ai^), at low energies 
as 

Vc,.A{r) = I Kc--iv(r - r')p(r')rfV, 



where p{r) is the nucleon distribution in the nucleus of nucleon number A and Vcc-N{f) the 
charmonium-nucleon potential. Since the density in the central core of nuclei is practically 
constant, falling sharply at the edges, and the range of Vcc-n{j) is very short compared to 
the size of nuclei, the above expression for Vcc-a{''^) is approximated by the nuclear matter 
result 

Ks-A = Po / V,-Mr)d'r = = -11 MeV, 
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for po = 0.17 fm ^ and m^ed the cc-N reduced mass, a result consistent with Ref. [jT2|, and 



comparable with the binding energy of protons and neutrons in nuclear matter. For finite 
nucleus, we describe the nuclear density by the Fermi density function p(r) = pQ/{l+e^^~'^^^'') 



which fits very well the data and incorporates the thickness of the nuclear surface [^. For 



the light nuclei we have used various forms for p(r) found in the literature to describe the 



data . We have investigated the binding of charmonium to various nuclei using different 
forms for Vcc-N{f) which reproduce the same scattering length = —0.2 fm for various 
potential ranges. Consistent results were obtained for Yukawa, Gaussian, and Bargmann 
potentials. As an example, we give in Table | the results obtained for a Gaussian form 
V{r) = l/oe-^"/^)', with a, = y?m,edK)^V2. 

We have determined the low energy interaction of heavy quarkonium and nucleons in a 
model independent way using the operator product expansion and the trace anomaly. The 
value of —0.2 fm for the cc-nucleon scattering length corresponds to a total cross section 
of about 5 mb at threshold. This value reproduces the anomalous behavior of the spin 
correlation A^n observed at ^/s = 5 GeV in pp elastic scattering at 6*0.111. = 90° near the 
charm threshold, by the enhancement ofpp amplitudes from the unitarity condition. The low 
energy charmonium-nucleon interaction only depends on the short range gluonic exchange 
properties of QCD between color singlet objects and do not involve meson degrees of freedom. 
Binding of charmonium was investigated for light, medium and heavy nuclei using realistic 
nucleon distributions for various forms and ranges of the cc-N potential. Our results show 
that nuclear bound states of charmonium could be formed for A > 6. A value of about 
10 MeV is found for the binding of charmonium in nuclear matter, which is comparable with 
the binding energy of nucleons, a remarkable result since the nucleon-nucleon interaction is 
two orders of magnitude stronger than the cc-N interaction. This is a consequence of the 
absence of Pauli blocking in the charmonium-nucleon system. Measurement of the J/ip- 
nucleon scattering near threshold is important to determine the van der Waals strength 
since this reaction is dominated by gluon exchange [0]. Study of charmonium production 
close to threshold p6| and its interaction with nuclear matter [|1^ would provide a better 
understanding of the mechanisms of color dynamics. 
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FIGURES 

FIG. 1. Ann as a function of piab at ^c.m. = 90° near the charm threshold, where the model is 
applicable. The data is taken from Ref. 4. The value of PQCD alone is 1/3. The upper, middle and 
lower curves correspond to = 0.6 fm, re = 0.8 fm, and re = 1.2 fm respectively and = —0.2 fm. 
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TABLES 



TABLE L Binding energies {H) of the r]c to various nuclei for a variational calculation cor- 



responding to realistic nuclear densities |25| and a Gaussian form for the cc-N interaction. The 
values of {H) are in MeV, the charmonium-nuclei reduced masses, Mj-ed) are in GeV and the range 
of the potential, R, is in fm. 





Mrcd 




(-H')_R=0.8 




^He 


L66 


> 


> 


> 


^Li 


L95 


-0.12 


-0.07 


> 


^Be 


2.21 


-1.31 


-1.04 


-0.68 


12c 


2.36 


-2.52 


-2.02 


-1.75 




2.44 


-3.31 


-2.92 


-2.54 


«Ca 


2.77 


-6.13 


-5.83 


-5.31 


56pg 


2.83 


-6.70 


-6.48 


-6.23 


208pb 


2.95 


-9.24 


-9.20 


-9.12 
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